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Abstract 

The perceived lateral position of a transmitted beam, upon propagating through a slab 
made of homogeneous, isotropic, dielectric material at an oblique angle, can be controlled 
through varying the velocity of the slab. In particular, by judiciously selecting the slab 
velocity, the transmitted beam can emerge from the slab with no lateral shift in position. 
Thereby, a degree of concealment can be achieved. This concealment is explored in numerical 
calculations based on a 2D Gaussian beam. 
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1 Introduction 

The topic of invisibility — which is a very old one in optics and electromagnetics (Wolf & 
Habashy 1993) — has lately acquired a new lease of life with the advent of metamaterials 
(Fedotov, Mladyonov, Prosvirnin & Zheludev 2005; Alii & Engheta 2006). In a similar vein, 
it has recently been proposed that the exotic electromagnetic possibilities offered by metama- 
terials may be harnessed to achieve cloaking, at least in principle (Leonhardt 2006; Milton & 
Nicorovici 2006; Pendry, Schurig & Smith 2006). The theoretical arguments underlying this pro- 
posed cloaking are based on the facilitation of coordinate transformations by nonhomogeneous 
metamaterials. 

A quite different approach to concealing a material is pursued in this paper. It is based 
on the perceived deflection of light by a material slab translating at constant velocity. No 
special material properties are required, but for illustrative purposes, we consider an isotropic, 
homogeneous, dielectric material. In two previous studies, we have demonstrated that the 
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much-heralded negative— phase-velocity phenomenon often associated with negatively refracting 
electromagnetic metamaterials can be realized by conventional materials through the process of 
uniform motion (Mackay & Lakhtakia 2004, 2006a). Here, we demonstrate that a substantial 
degree of concealment may also be realized by uniform motion. 

As regards notational matters, 3 vectors are in boldface, with the " symbol identifying unit 
vectors. Double underlining signifies a 3x3 dyadic and / is the identity 3x3 dyadic. The 
superscript denotes the transpose of a column vector. The permittivity and permeability of 
vacuum are Cq and fiQ. The vacuum wavenumber is ko = LO-^/e^JI^ with cj being the angular 
frequency, and the vacuum wavelength is Aq = 2'K/ko. 



2 Planewave propagation into a uniformly moving half— space 



As a preliminary to concealment of a moving slab (§3), let us consider a uniformly moving half- 
space. Suppose that a plane wave is launched with wavevector kj = fcjkj from vacuum {z < 0) 
towards the half-space z > occupied by an isotropic, nondissipative, dielectric material. This 
material moves at constant velocity v = f v = vx, parallel to the interface and in the plane of 
incidence. In an inertial frame of reference that moves with the same velocity v with respect to 
the laboratory frame of reference wherein kj is specified, the refracting material is characterized 
by relative permittivity e^. The Minkowski constitutive relations of the moving half-space in 
the laboratory frame of reference are (Chen 1983) 



T>{r) = eoer^'E{r) + ^/^{mx I) • H(r) 
B(r) = -,/e^(m X I) • E(r) + //q^ • H(r) 



(1) 



where 



a = a / + (1 — a) vv , 

/3 (e, - 1) 



a = 



m = mv , m 



,2 , l3 = Vy/€^a 



> . 



(2) 



1 - 

In order to exclude the possibility of evanescent plane waves, > 1 is assumed. The envisaged 
scenario is illustrated schematically in Figure 1. 

The angle between the refracted wavevector kj = fetk^, as observed from the laboratory 
frame of reference, and the unit vector z normal to the interface is related to the angle of 
incidence Ai = cos' -1 (ki • z) by (Chen 1983) 



bt = sm 



I f ko sin ( 



where 



kf. — kn 



;o|l + c[l-/3(ki -v) 



1/2 



(3) 



(4) 



is the wavenumber of the refracted wave and ^ = (e,- — 1)/(1 — P"^)- Since < < 7r/2V( 
(0,7r/2), refraction is positive G (-1, 1) (Mackay & Lakhtakia 2006b). 
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The time-averaged Poynting vector of the refracted plane wave is given by (Chen 1983) 

Pt = Pt = ( |Ci|2 + e,|C2p) (kt X v)2 [kt + CP (K - pkt • v) v] , (5) 
where Ci and C2 are constants, and the angle between z and Pt is 

cl>p = tan-i ( I . (6) 

As an illustrative example, the angle <^p is plotted in Figure 2 against /? G (—1,1) for 0, G 
{15°, 45°, 75°}, for the half-space characterized by = 6.0. The orientation of the refracted 
time-avcragcd Poynting vector rotates towards the direction of motion as (3 increases from — 1. 
The counterposition regime — which occurs where 0p < for 0t > — is discussed elsewhere 
(Lakhtakia & McCall 2004; Mackay & Lakhtakia 2006b). 

In connection with Figure 2, it is of particular interest here that 0p = 4>i at (i) /? = 0.08 for 
(t)i = 15°, (ii) P = 0.29 for = 45°, and (iii) P = 0.78 for = 75°. That is, there exist angles of 
incidence at which the time-averaged Poynting vector is not deflected by the uniformly moving 
half-space. This suggests that it may be possible for a light beam — not to be confused with a 
plane wave — to pass through a uniformly moving slab at an oblique angle without experiencing 
a lateral shift in position. That suggestion inspired the research presented in the next section. 



3 Beam propagation through a uniformly moving slab 



Suppose that the uniformly moving half-space considered in §2 is now replaced by a slab of 
thickness L moving at constant velocity v = vx parallel to its two surfaces, as schematically 
illustrated in Figure 3. The slab — which is characterized, as before, by > 1 in a co— moving 
reference frame — is immersed in vacuum. 



A 2D beam with electric field phasor (Haus 1984) 

/oo 
ei(i?)*(i?) exp[i(ki - r)] di9, 
-00 



z<Q 



(7) 



is incident upon the slab at a mean angle 9i relative to the slab normal direction z. The beam 
is represented as an angular spectrum of plane waves, with 



kg 



(^^ cos 9i + Vl-^^ sin Oi^ x - (^t? sin Oi - Vl - cos 



(8) 



being the wavevector of each planewave contributor. The angular-spectral function ^'(i?) is 
taken to have the Gaussian form (Haus 1984) 



w(w) = ^_ exp 



2tt 



{ko Wo -dy 



(9) 



with Wq being the width of the beam waist. Two polarization states are considered: parallel to 
the plane of incidence, i.e., 



ei{-&) = ey = (^^ smOi - ^1-^^ cos^^) i + cos^j + ^1-^^ sin^j) 



(10) 
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and perpendicular to the plane of incidence, i.e., 



ei{i}) = e_L = y. 



(11) 



The electric field phasor of the reflected beam is given as 

/oo 
er{'&)'^{'&) exp[i(k^ • r)] d^, 
-oo 



z < 0, 



with 



and 



1? cos 9i + sin 6i \ X -\- ("d sin cost/ 



- (^i? sin 0i - Vl-i?2 cos ic 
+ (^?? cos 0, + Vl-??2 sin z 



r_L e_L 



for ei{'&) = ey 
for ei{-d) = e_L 



(12) 
(13) 

(14) 



The electric field phasor of the transmitted beam is given as 

/oo 
et(??)*(??) exp{i[kt • (r-Lz)]} di9, 
-oo 



z > L, 



(15) 



with kf = kj and 



t|| ej(i?) for ej('!?) = ey 
for ej(i?) = e_|_ 



Expressions for the reflection coefficients r||^_L and transmission coefficients iy 
equations (34)-(37) in the Appendix. 



(16) 

are provided in 



In view of Figure 2, we fixed the mean angle of incidence of the beam at 9i = 45° and 
explored the behaviour of the transmitted beam for /? < 0.29, /? = 0.29 and /? > 0.29. The 
energy density of the beam in both half-spaces, as measured by 



lEl 



|Ei + EJ2 for z<0 



IE. 



for z > L 



(17) 



is mapped for z/A,, G (—8, 12) and .t/A„ G (—25, 25) in Figure 4 with the slab thickness L = 4Ao. 
The restriction i? G [—1, 1] was imposed to exclude evanescence. A beam waist of Wq = 1.75Ao 
was selected for all calculations. We considered (3 G {—0.15,0.29,0.8} for both = ey and 
Bi = e±. 

Regardless of the polarization state, the transmitted beam does not undergo a lateral shift 
(relative to the incident beam) when /3 = 0.29. However, the transmitted beam is laterally 
shifted in the direction of — x when (5 < 0.29 and in the direction of +x when (5 > 0.29. The 
energy densities of the reflected and transmitted beams are sensitive to /3 and the polarization 
state of the incident beam. 

A more quantitative representation of the transmitted beam is provided in Figure 5, wherein 
|Ep is plotted against x f or /3 G {—0.15, 0.29, 0.8} at ^; = 4Ao. For comparison, |Ep for the beam 
in the absence of the moving slab is also plotted. It is clear that the beam position for /3 = 0.29 
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coincides with the beam position in the absence of the moving slab. At /3 = 0.29, the peak 
energy density of the transmitted beam for the case of parallel polarization is approximately 
11% less than it would be if the slab were absent; the corresponding figure for perpendicular 
polarization is 38%. 

The median shift of transmitted beam in relation to the incident beam is defined as 
A = (^J x\Et{x,L)\^ dx^ (^J \Et{x,L)\^ dx^ 

-(^j x\Bi{x,L)f dx^ (^J \Bi{x,L)fdx^ . (18) 

For both parallel and perpendicular polarizations, A is plotted against /? G (— 1, 1) in Figure 6. 
Thus, regardless of the polarization state, the beam can be shifted laterally along itx by means 
of uniform motion. In particular, the zero beam shift at /3 = 0.29 is further confirmed in 
Figure 6. 



4 Concluding remarks 

Our numerical investigations show that a 2D beam can pass obliquely through a uniformly 
moving slab without undergoing a lateral shift in its position. At a fixed angle of beam incidence, 
this effect occurs only for a unique translational slab velocity. However, extrapolating from 
Figure 2, for every angle of beam incidence a slab velocity can be found at which the beam 
undergoes no lateral shift. Furthermore, for a fixed angle of incidence, a pulsed beam will 
undergo zero lateral deflection, provided that the constitutive parameters do not vary with 
angular frequency in the pulse spectrum. 

The degree of concealment achieved by uniform motion is not 100% due to reflections but, 
in the particular case of the example considered in §3, almost 90% of the peak energy density 
of the beam can be transmitted without deflection. 
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Appendix 

The reflection coefficients ry _|_ and transmission coefficients are straightforwardly calculated 
by solving the reflection-transmission problem as a boundary value problem. We outline the 
procedure here, further details being available elsewhere (Lakhtakia Sz Messier 2005). 

Consider the plane wave with electric and magnetic field phasors 

Ei{x,z) = e{z,6) exp{ikoXsin9) 
H{x,z) = h.{z,6) exp (zfcoX sin 0) 
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propagating in the xz plane. As in §3, a moving slab described by the Minkowski constitutive 
relations (1) occupies the region between z = and z = L; elsewhere there is vacuum. We write 

piz,e)=p^{z,e)yi + pyiz,e)y + p,iz,e)z, {p = e,h). (20) 

Substitution of equations (1) and (19) into the source-free Maxwell curl postulates 



V X E(x, z) — iui'B{x, z) = 0, 

V X H(a;, z) + iujT>{x, z) = 0, 



(21) 
(22) 



delivers four differential equations and two algebraic equations. The latter two equations are 
easily solved for and hz- Thereby, the four differential equations may be expressed in matrix 
form as 

^ [i{z,d)] = iK[v{e)mz,e)], 



where 



is a column vector and 



dz 



i{z, 9)] = e^iz, 9), ey{z, 9), hx{z, 9), hy{z, 9) 



P(^) 



r?oP 

-770 

-e,.p/r/„ 

erho 



is a 4x4 matrix with 



a 



(m + sin^)^ 



The solution to (23) is conveniently expressed as 

[f(L,^)] = [M(L,e)] [f(0,^)], 

in terms of the transfer matrix 

[m{L,9)]=e^v{iK [P(^)] L}, 



(23) 
(24) 

(25) 

(26) 
(27) 
(28) 



Now we turn to the incident, reflected and transmitted plane waves. Let the incident plane 
wave be represented in terms of linear polarization components as 



ei[z,9) = [a_L y + a|| (— cos X + sin6' z)] exp (ifco^ cos 0) | 
hj(z, 6) = r]^^ [aj_ (— cos^x + sin^z) — ay y] exp (ikQZC0s9) J 

The corresponding reflected and transmitted plane waves are given as 

er{z,9) = [aj_r_L y + a||r|| (cos x + sin z)] exp (— ifco^; cos 0) | 
hr(2, 6) = rjQ^ [a±r± (cos x + sin6' z) — a||r|| y ] exp {—ikoZ cos 9) J 

and 

et{z, 9) = [a±t± y + a||t|| (— cos^x + sin^z)] exp [iko{z — L) cos 9] 
h.t{z, 9) = rjg^ [a±t± (— cos^x + sin^z) — a||t|| y ] exp [iko{z — L) cos9] 



z < 0. 



(29) 



z < 



(30) 



z>L, (31) 
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respectively. By application of the boundary conditions at z = and z = L to the solution 
(27), the reflection and transmission coefficients are found to be related by the matrix algebraic 
equation 

(32) 



[K(0)] [t^,t||,0,0]^ = [M{L,e)] [K{9)] [l,l,r^,r^^f, 



wherein 



Kie) 





1 

-tJq^ cos ( 




— cos^ 






1 

rjg^ cos ( 



cos 6 







Thus, after some manipulation, the reflection and transmission coefficients emerge as 



cos 



m = 



erp) sin (kf,L^erp) 



t± 



(cos^ 9 + Erp) sin (/c,)Ly^e^) + 2i^erp cos (kaLyJerp) cos 9 ' 

(p — e,f cos^ 0) sin {k^L^erp) 
(e^ cos^ 9 + p) sin {k^L^erp) + 2i^fe^ cos {k^L^t^p^ cos ' 
—2i.^/e^ cos 9 

(cos^ + e^p) sin (^kgL^erp) + 2i,^/ePp cos (fcoLy^e^) cos ' 

2iy^e,.p cos 

(e^ cos^ 9 + p) sin {koLy/€rp) + 2i^ej.p cos {k^L^erp) cos ^ 



(33) 



(34) 
(35) 
(36) 
(37) 
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Figure 1: A plane wave with wavevector kj is incident from vacuum onto a half-space occupied 
by a simply moving dielectric material at an angle (pi with respect to the unit vector z normal 
to the planar interface. The moving material is characterized by relative permittivity e,. > 
in a co-moving frame of reference. As observed in the non-co-moving (laboratory) frame of 
reference wherein the incident plane wave is specified, the refracted wavevector makes an 
angle (pt with z. 




Figure 2: The angle (j)p (in degree) between the time-averaged Poynting vector and the unit 
vector z, plotted as a function of (3 £ (—1, 1), when the angle of incidence (pi = 15° (solid curve), 
45° (dashed curve) and 75° (broken dashed curve); and = 6.0. The red lines indicate where 
(pP = 4>i- The counterposition regime {(pp < 0°,(pt > 0°} is shaded. 
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Figure 3: A beam is incident onto a simply moving slab at an angle 9i with respect to the 
unit vector z normal to the planar interface. The moving material is characterized by relative 
permittivity > in a comoving frame of reference. As observed in the non-comoving (labo- 
ratory) frame of reference wherein the incident plane wave is specified, the transmitted beam is 
shifted by A, parallel to x, relative to its position if the slab were absent. 
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Figure 4: Normalized |Ep is mapped in the xz plane for a 2D Gaussian beam incident onto 
a simply moving slab at an angle 9i = 45°. The relative speed of the slab is: /3 = —0.15 
(top); P = 0.29 (middle); and (3 = 0.8 (bottom). The electric field phasor of the incident beam 
is polarized parallel (left) and perpendicular (right) to the plane of incidence. The red line 
indicates the mean beam position in the absence of the moving slab. 
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Figure 5: Normalized |Ep at z = 4Ao for /? = —0.15 (broken dashed curve); /? = 0.29 (solid 
dark curve); and /? = 0.8 (dashed curve). The solid red curve represents the normalized |Ep 
in the absence of the moving slab. The electric field phasor of the incident beam is polarized 
parallel (left) and perpendicular (right) relative to the plane of incidence. 




Figure 6: The median beam shift A at z = 4Ao, plotted against /3. The electric field phasor of 
the incident beam is polarized parallel (left) and perpendicular (right) relative to the plane of 
incidence. Zero median beam shift at /? = 0.29 is indicated by dashed lines. 
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